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Abstract. The entanglement of two qubits, each defined as an effective two- 
level, spin 1/2 system, is investigated for the case that the qubits interact via 
a Heisenberg XY interaction and are subject to decoherence due to population 
rclaxation and thermal effects. For zero temperature, the time dependent 
concurrence is studied analytically and numerically for some typical initial states, 
including a separable (unentangled) initial state. An analytical formula for non- 
zero steady state concurrence is found for any initial state, and optimal parameter 
vàlues for maximizing steady state concurrence are given. The steady state 
concurrence is found analytically to remain non-zero for low, finite temperatures. 
We also identify the contributions of global and local coherence to the steady 
state entanglement. 



PACS numbers: 03.65.Yz, 03.67.Mn, 75.10.Jm, 05.50.+q, 42.50.Lc 

1. Introduction 

Entanglement is a property of corrclatcd quantum systems that cannot be accounted 
for classically. Entangled states of distinct (possibly interacting) quantum systems, 
which are those that cannot be factorized into product states of the subsystems, 
are of fundamental interest in quantum mechanics. The production of pairwise 
entangled states is an essential requirement for the operation of the quantum gates 
that make quantum information and quantum computation possible £Q. Considerable 
attention has been devoted to interacting Heisenberg spin systems [21 El EU , 
which serve as a model for various sòlid state □ E] or NMR [TÜ] quantum 
computation schemes and for simulating magnètic phenomena in condensed matter 
systems using àtoms in optical lattices |1 11 112) . Indeed, general Hamiltonians that 
include Heisenberg spin-spin interactions have been proposed as "genèric" ^SJ or 
"ideal" ^1] model Hamiltonians for quantum computation systems. A key question 
for entangled quantum states is the effect of decoherence due to the environment (see, 
e.g., |15l 1161 1171 1181 [TÜ] and references therein), which is not only a fundamental 
issue for quantum computation devices |2()l l?Tj but also for the relation between 
quantum and classical physics ^3 • Although there have been many investigations 
of decoherence in recent years, careful investigation of well-understood model systems 
continue to produce surprises that add to fundamental understanding. For example, 
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Yu and Eberly |2.'ij have recently shown that the entanglement of a pair of non- 
interacting qubits in the presence of spontaneous decay of the upper states may 
decohere in a finite time instead of exponentially. 

In this paper we examine decoherence due to both population relaxation and 
thermal effects for an entangled (and interacting) two qubit system. The Hamiltonian 
for our two-qubit system has the form of the well-known Heisenberg XY model for 
two interacting spins in the presence of an external magnètic field, where the effective 
magnètic field is defined by the energy separation of the two-level system that we 
associate with each (spin 1/2) qubit. As noted above, this form of Hamiltonian is very 
common in models for quantum computing |HM12llTÏÏlll4| . Our analysis of decoherence 
complements that of Ref. (22| by examining a system in which the qubits interact. For 
our two-qubit model system at zero temperature, we find that for any initial state, 
including the common one in which the two qubits are initially unentangled, the system 
reaches a steady state of pairwise entanglement in spite of population relaxation. 
The extent of steady-state entanglement is sensitive to both the spatial anisotropy of 
the interaction between qubits and to the energy level separation of the two levels 
associated with each qubit. To the extent that these two parameters can be varied in 
some particular physical realization of our model system, the magnitude of steady state 
entanglement may thus be controlled. We analyze both analytically and numerically 
the time-dependent evolution of the entanglement (as measured quantitatively by the 
concurrence of our model two-qubit system for some typical initial states: a 

pure, separable initial state; a pure, entangled initial state; and a mixed initial state. 
We also obtain an analytic formula for the steady state concurrence that shows its 
dependence on both the system parameters and the decoherence rate and that enables 
us to specify optimal vàlues for these parameters to achieve the maximum possible 
concurrence. In a separate section, we consider the case of finite temperature and 
present an analytic formula for the concurrence, which remains non-zero over a finite 
range of low temperatures. In our concluding section, we discuss some implications of 
these results. 

2. Two-Qubit Hamiltonian 

We note that the Hamiltonian of a Heisenberg chain of N spin i partides with nearest- 
neighbor interactions is 



where S" = ^c™(a = x, y, z) are the local spin i operators at site n, the er™ operators 
are the Pauli matrices at site n, the periòdic boundary condition Sm+i = S\ applies, 
and h = 1. For arbitrary J Q 's, the Heisenberg chain is often called the XYZ model. 
The chain is said to be antiferromagnetic for J a > and ferromagnetic for J a < 0. 
The XY(J Z — 0) and the Heisenberg-Ising (J y = J z = 0) interactions have been 
analyzed for nuclear spin systems |Ü], in particular for nuclear magnètic resonance 
approaches to quantum computation (see, e.g., Section 7.7 of Ref. |10|1. 

The Hamiltonian H for an anisotropic two-qubit Heisenberg XY system in an 
(effective) external magnètic field u> along the z-axis is: 



JV 




(D 



H = üj{SI + 5 2 2 ) + J{S+S^ + 5f S+) + A{S+S+ + S^S^) 



(2) 



Entanglement Evolution in the Presence of Decoherence 



3 



where J = (J x + J y )/2, A = (J x — J y )/2, and S ± — S x ±iS v are the spin raising 
and lowering operators. The first term in the Hamiltonian describes the energy of 
the spins in the efFective external magnètic field. This effective field is dehned by 
the energy levels of our qubits: we assume that each of our two qubits represents 
an identical two level system whose two energies are defined by ±lu/2. The spin 
interaction Hamiltonian, described by the second and the third terms in Eq. P|. 
produces the coherence of the two qubits that is necessary for their entanglement in 
the presence of decoherence. As shown below, the third term, whose magnitude is 
proportional to the parameter A, which describes the spatial anisotropy of the spin- 
spin interaction, is essential for the production of steady state entanglement. 

3. Time Evolution of the Concurrence at Zero Temperature 

The time evolution of the system for zero temperature, T — 0, is given by the following 
master equation (see, e.g., [2ül|27] and Section 8.4.1 of |l(Jjh 

p = -i[H,p]+1D[Sï]p + 1D[Sï]p. (3) 

Here p is the density matrix, which in the presence of population relaxation represents 
the mixed state of the system. The Lindblad súper operator T> '28 k is defined by 
T>[A]B = ABA^ — {A^ A, B} /2, which describes the population relaxation of the 
upper state of each qubit due to the environment; 7 is the phenomenological rate of 
population relaxation, which for simplicity is assumed to be the same for each of the 
two qubits (i.e., we assume each qubit has the same interaction with the environment). 
As discussed below, the assumption of a single decay rate, 7, requires us to place 
restrictions on the magnitude of the coupling between qubits. 

Entanglement is incrcasingly regarded as a physical resource of a quantum 
information system (see, e.g., Section 12.5 of Ref. |1U)1 and many measures for 
quantifying entanglement have been developed (see, e.g., [24, 25, 29 l30ll3*Tll32p . Since 
decoherence processes cause the system state to become mixed, we use the measure 
of entanglement termed concurrence |24l 125) . For a system described by the density 
matrix p, the concurrence C is 

C = max(v/XI"- V^2 - \J% - V^O) , (4) 

where Ai, A2, A3, and A4 are the eigenvalues (with Ai the largest one) of the "spin- 
flipped" density operator R, which is defined by 

R = p {(7y ® <Ty) p* (a y ® (Xy) , (5) 

where p* denotes the complex conjugate of p and a y is the usual Pauli matrix. C 
ranges in magnitude from to 1; nonzero C denotes an entangled state. 

The basis states \ij>i) for our two-qubit system are the separable product states 
of the individual qubits: 

|V>i> = |e)i® \e) 2 , 

\ih) = |e)i <8> 15)2, 

1^3) = \g)i® |e} 2 , 

1^4) = \g)i ® \g)i. (6) 
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In general, a two-qubit system is represented by a density matrix having sixteen non- 
zero elements. For our Hamiltonian, however, the density matrix can be represented 
as the sum of two submatrices that evolve independently of one another, 



P = 



( Pn 



V P41 





P22 
P32 






P23 
P33 
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P21 



P31 




Pl2 




P42 



Pl3 
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\ 



(7) 



Pu \ 



Pi4 ) 

i.e., in solving Eq. for p(t) the forms of each of the two submatrices in Eq. @ are 
preserved. (Note that the second matrix on the right hand side of Eq. Q does not 
have the form of a density matrix.) Each of the cxamples in this paper has an initial 
density matrix whose form is that of the first matrix on the right of Eq. Q). This 
limitation is not very restrictive, as unentangled, entangled, and mixed states can all 
be described. Furthcrmorc, for a state having a density matrix of the form of the first 
matrix on the right of Eq. (jJJ), the concurrence has the following analytic form: 

C = max{0,Ci,C 2 }, (8) 

where 



Ci = 2(|p4i| - VP33P22) 

C2 = 2(|p 3 2|-vpïïpir)· (9) 

The solutions of the master equation in Eq. simplify by transforming from the 
product state basis \ipi) in Eq. © to the basis of eigenstates of the Hamiltonian 
in Eq. 101, 

A 



1*2) 
1*3) 
|$ 4 ) 

n 

7V± 



^V + (l.95> 
1 

1 

N-(\gg) 



feg) 
-À\ge) 



Q — uj 

\ge)), 
\eg}), 

A 

n + uj 



\ee)), 



\ee)), 



Vu 2 + A 2 , 

(OT W )/VA 2 + (fiTw) 2 . 



(10) 

(11) 
(12) 



After transformation, the solutions for each element p aa > of the density matrix (where 
p denotes p in the eigenstate basis) can be found analytically. For the interesting 
special case that both qubits are initially in their ground states (i.e., the system is 
initially in state ^4) in Eq. (JüJ), the analytic solutions for p aa '(t) are: 



1 



— uja + 



2Üa 

+ íl{a - 2A 2 ) + 2e^ 7t A 2 7 sin [20í] 
A 2 r 

P22W) = 7t— O — ne- 2jt - e- 7í 7 sin [20íl 
Sia - 

P33{t) = p2l{t) 

pa{t) = 1 - p U (t) - p 2 2(t) - P33(t) 



(13) 

(14) 

(15) 
(16) 
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; e - 7 * cos [2ííí] 



+ 20e~ 7 * sin [20í] + 7 



(17) 
(18) 



P4l(t) = Puí*) 



where all other matrix elements are zero and where 



a = 4Í7 2 + 7 2 . 



(19) 



From Eqs. (|13I18|) it is evident that both the off-diagonal (coherence) matrix 
elements pu and pn in Eqs. (|17I18|) and the diagonal (population) matrix clements 
p aa in Eqs. (|13I16|) have terms that oscillate with frequency 20. Note that 
the coherence matrix elements pn and pn are non-zero only when the spin-spin 
interactions are anisotropic (i.e., A ^ 0); also, the value of fí is sensitive to this 
anisotropy (cf. Eq. From Eqs. (|13I18|) it can be seen that the coherence matrix 

elements have terms that decay at the rate 7 while the population matrix elements 
also have terms that decay at the rate 27. Analytic solutions similar to Eqs. I|13ll8fl 
can be given for some other initial states. 

The assumption of a single decay rate, 7, in the master equation (0 requires some 
discussion. Owing to the interaction between qubits described by the Hamiltonian 
J2Jl, the two-qubit energy level structure is altered from that describing non- 
interacting, identical qubits. Nevertheless, the assumption of a single decay rate, 
7, is reasonable provided the interaction does not significantly alter the effective 
energy level separations, or, more precisely, provided the rotating wave approximation 
remains vàlid [SS] (see, e.g., pp. 160-161 of Ref. The eigenstates in Eq. (fTTTft have 

the following eigenenergies 0]: the Bell states have eigenvalues ±J while the other 
eigenstates have eigenvalues ±Í2. Thus if we restrict the magnitudes of the coupling 
parameter J and the anisotropy parameter A to vàlues such that, 



we shall ensure that the energy level separations of the interacting qubit system do 
not change by more than 10% from that of the non-interacting qubit system. Except 
where it is explicitly mentioned otherwisc, all cxamples given below have parameter 
vàlues for which the above inequalities are satisfied. 

Perhaps surprisingly, the decoherence due to population relaxation does not 
prevent the creation of a steady state level of entanglement, regardless of the initial 
state of the system. This is demonstrated in Figs. ^ and [21 which show the time 
evolution of the concurrence (cf. Eq. (0}) for three different initial states: (1) An 
unentangled, separable state, (cf. Eq. ijïïjl); (2) a completely entangled state, the 
Bell state \ty) = -^(\gg) — |ee)); and (3) a mixed state, defined as an equal mixture of 
\ipi) and the Bell state | < E > 2)· In Fig.^we consider the case that J — A = w/10, which 
implies that J y = and which thus corresponds to the "genèric" quantum computation 
model Hamiltonian of Ref. In Fig.|21we consider the case that J = w/10 and that 
A = 0.458a;, which corresponds to a general case in which J x and J y have oppositc 
signs, which may possibly be achieved for an optical lattice system For each 

of the three initial states considered, the corresponding curves in Figs. ^ and [3 give 
the numerical results for the concurrence defined by Eq. (@J, after solving Eq. © 
numerically for the density matrix in the separable representation (cf. Eq. I©)- Since 
each initial state has a density matrix of the form of the first matrix on the right of 



|J|/w < 0.1 
{ü-üj)/üj < 0.1, 



(20) 
(21) 
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Figure 1. Plot of T = concurrence vs. scaled time, yt, for three different initial 
states: (1) An initially unentangled state, |\l>) = \gg) (sòlid line); an initially 
entangled state, the Bell state |<P) = -^(\gg) — |ee)) (dashed line); and (3) an 

initially mixed state, defined as an equal mixture of \gg) and -j= (\eg) + \ge)) (sòlid 

squares). The system parameters are: 7 = 0.3, lo = 1.0, J = 0.1, and A = 0.1. 



Eq. (0, the concurrence for each of these states is given also by Eqs. í|819|> . (Note that 
discontinuities in the time derivatives of C(t) for the dashed curve in Fig. in the 
range 2.0 < t < 2.5 stem from the dcfinition in Eq. (0}; all dcnsity matrix elements 
are smooth functions of t.) The sòlid circles on the curves for the initial state \gg) 
in Figs. n and give the concurrence obtained from the analytic Eqs. (|8I9J) (after 
transforming the analytic expressions in Eqs. - (|18|l for this state's density matrix 
paa' to pij). These analytic results coincide with those obtained by direct numerical 
solution of Eq. 10). 

Despite the presence of decoherence, the results in Figs. ^ an d [21 show that the 
concurrence reaches the same steady state value (after some oscillatory behavior) 
for a given set of system parameters, regardless of the initial state of the system. 
(This is true even for initial states having non-zero matrix elements belonging to the 
second matrix on the right of Eq. J7J; for our system, such matrix elements vanish 
in the steady state.) Clearly the Heisenberg spin-spin interaction in Eq. J5J) serves to 
maintain an entangled state despite the presence of decoherence in Eq. J2J. We find a 
steady state concurrence of 0.09309 for the system parameter vàlues chosen in Fig. ^ 
and a steady state concurrence of 0.28916 for the system parameter vàlues chosen in 
Fig. |3 

The analytic expressions for the T = steady state vàlues of p%j(t) are as follows: 

A 2 

Pll = P22 = P33 = — (22) 
Q 

3A 2 

p 44 = 1 (23) 
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Figure 2. Plot of T = concurrence vs. scaled time, ■yt, for the same three initial 
states as in Fig.0 but for the following differcnt system parameters: 7 = 0.458, 
ui = 1.0, J = 0.1, and A = 0.458. 

-2cjA-iA 7 

Pia = (24) 

a 

Pa = Pia ( 25 ) 
The corresponding T — steady state concurrence is found to be: 



_ 2VA2(4^ +7 2)-2A 2 

Csteady — — (,^0 ) 

This result stems from C\ in Eq. JÜJ calculated for the separable basis density matrix 
Pij in Eqs. I|22I25(1 . The steady-state concurrence is seen to depend on the system 
parameters u>, A, and 7 (but not on J). Also, 7 serves as a scale factor, i.e., C s teady 
depends only on the scaled variables, Q — w/7 and A = A/7. These parameters 
may be varied in order to maximize Csteady- The function Csteady {ü>, A) is shown in 
Fig. |21 one sees that the surface has a ridge along which it takes its maximum value. 
The coordinates of the ridge and the value of C s tead y on the ridge may be determined 
analytically. For fixed u>, C s teady (cf. Eq. (|26fl ) has its maximum at the following 
value of A: 



V4^ + 7 2 

max ~ (1 + V5) ' [ } 

The sòlid line in the ü-A plane of Fig. [21 represents the locus of points A max (üi) 
given by Eq. (|27|l (upon division by 7). Substitution of A max into Eq. l|26|l gives the 
parameter-independent maximum value of the concurrence, represented by the sòlid 
line in Fig. along the ridge of Csteady'- 

Csteady{A max ) = (1 + Vs)" 1 = 0.309. (28) 

Eq. I|2t)|) shows that in order to have a positive value of Csteady-, one must have 
4uj 2 + 7 2 > A 2 . Note finally that Eq. (j2"5jl was derived from Eq. ltÏÏ5|) without taking 
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Figure 3. Plot of the T = steady state concurrence (cf. Eq. 1261 ) as a function 
of the scaled energy ü and the scaled anisotropy parameter À (ranging from 0.309 
to 1), where (D = uj/·y and À = A/7. The sòlid lincs locate the maximum value 
of concurrence (cf. Eq. 1281 *): sec text for discussion. 



into account the restrictions on the parameter vàlues imposed by the conditions in 
Eqs. Q20|l and (|21[) that are necessitated by our assumption of a single decay rate, 7. 
Nevertheless, one sees for the example plotted in Fig.[3that there do exist vàlues of 
the parameters that satisfy Eqs. I|2U[1 and i|21|) for which one obtains a steady state 
level of concurrence that is close to the global maximum value given by Eq. (|28|) (and 
shown by the sòlid line in Fig.[2|). 

4. Temperature Dependence of the Steady State Concurrence 

It is of interest to examine how the steady state entanglement obtained for zero 
temperature in the prior section changes when the temperature is finite. For simplicity, 
we assume that each qubit interacts with the same thermal bath. It is known that 
the equilibrium entanglement must vanish at some finite temperature [34] . In order 
to examine the effect of thermal decoherence on the entanglement for our system we 
consider the following master equation EÏÏ>| : 

p = - i [H, p] + 7 (n + 1)23 [Sf] P + 7Í« + ^ [5 2 ~] p 

+ jnV [S+] p + jnV [S+] p, (29) 

where fi, the average excitation of the thermal bath, parametrizes the temperature. 
Note that n is zero at zero temperature, whereupon one observes that Eq. (|29|l reduces 
to Eq. 10; also, n becomes infmite as the temperature becomes infinite. The master 
equation l|29|) may be solved to obtain the following analytic expressions for the steady 
state density matrix of our system: 

_ n 2 (4ïü 2 + (1 + 2nf ) + A 2 (l + 2n) 2 
911 ~ (1 + 2n) 2 (4d> 2 + (1 + 2n) 2 + 4A 2 ) 
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1 M 4w 2 + (1 + 2n) 

P22 = P33 = t[1 - 



P44 



4 L (l + 2n) 2 (4cj 2 + (l + 2n) 2 + 4A 2 ) J 
4w 2 (l + nf + (1 + 2n) 2 ((l + n) 2 + A 2 ) 
(1 + 2rl) 2 (4w 2 + (1 + 2n) 2 + 4A 2 ) 



À(2tD + i(2n+l)) 
Pl4 (l + 2n)(4<I> 2 + (l + 2n) 2 +4A 2 ) 1 ' 

In the limit of zero temperature (i.e., n — ► 0), the density matrix elements in Eq. (|30|) 
reduce to the results in Eqs. (|22|l - l|25|l . In the limit of infinite temperature (i.e., 
n — > 00), the density matrix becomes diagonal, with each diagonal element equal to 
1/4, indicating, as expected that all entanglement vanishes. 

The concurrence may be calculated for the finite temperature, steady-state 
density matrix in Eq. (|30[1 to obtain: 



v/A 2 (4<I> 2 + (1 + 2fï) 2 ) 1 



C(Q, A, n) - 2- (i + 2 _ )(4 - 2 + (i + 2 _ )2) 2 
(4w 2 + (l + 2n) 2 ) 



2[(l + 2n) 2 (4í] 2 + (l + 2fï) 2 )]' 

where all system parameters have been normalized by the relaxation rate 7: A = A/7, 
üj — w/7, and n = íl/7 = Vil" 2 + A 2 (cf. Eq. JHJ). In the limit of zero temperature 
(i.e., tt. — > 0) , the concurrence in Eq. fjHip reduces to that in Eq. (|26J) . The behavior of 
this finite temperature, steady state concurrence is shown in Fig.0]for the same two 
sets of system parameters considered in Figs.^and|21respectively. One sees that both 
curves decrease with increasing n until eventually the concurrence vanishes at a finite 
value of n, as expected jHJ- One sees also that the larger the value of the interaction 
asymmetry parameter A, the larger the value of the concurrence at any finite value of 
n. For any fixed temperature (i.e., n), as the effective magnètic ficld, üj, increases, the 
concurrence takes a finite, non-zero value. In the limit u> — * 00, one has that fi — * 
and C — > |A|/(D. This decrease with Q^ 1 as well as the fact that C > only if A 7^ 
is consistent with the results of Ref . 0] . 



5. Discussion and Conclusions 



Quantum coherence is a necessary requirement for the existence of entanglement. One 
may define coherence existing in a singlc qubit as "local coherence" while coherence 
between two qubits A and B may be defined as "global coherence" ■ How do local 
and global coherence relate to entanglement? The answer for our model system may 
be understood by considering the relation between a general two-qubit density matrix, 
Pab, and the reduced density matrices, p A and p B , for each of the two qubits, where 
p A = trB(p AB ) is obtained by tracing over the degrees of freedom of qubit B, and 
similarly for p B . In our model, after doing partial traces, we find that in the steady 
state the local coherence of each qubit is zero, i.e., p A and p B are diagonal matrices. 
However, there exist global coherence terms for our system (i.e., pu and pn) that are 
non-zero, indicating that global, not local, coherence is responsible for this system's 
entanglement. For time t > 0, pu and pn for our system are always non-zero; P23 
and P32 may be non-zero for finite times, but vanish in the steady state. Ref. 23 
considered a decohering system of two entangled (but non-interacting) qubits. In 
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Figure 4. Plots of the finitc tcmpcrature, steady state concurrence (cf. Eq. 1311 ) 
as a function of the average thermal excitation function, n, for the same two sets 
of system parameters as in Figs.0and[5] The curve in the lower left còrner of the 
figure corresponds to the same system parameters as in Fig. the curve close to 
the diagonal corresponds to the same system parameters as in Fig. Note that 
at n = both curves begin at the steady state vàlues of the concurrence shown 
in Figs.QandEl 



that system, both local and global coherence vanished in the asymptotic time limit; 
however, in some cases, the latter vanished for finite times |23| . 

An interesting question regarding the steady state of our model system is whether 
or not it is a decoherence- free subspace [37]. Typically, a decoherence- free subspace 
is defined to be one for which the decoherence terms in the system's master equation 
vanish |38| . For our system, this would mean that the second and third terms on 
the right of Eq. J2J vanish for the case of zero temperature or, for the case of finitc 
temperature, that all terms except for the first one on the right of Eq. I|29[l vanish. 
However, in our system, the decoherence terms in Eqs. (0 and (|29[1 do not vanish; 
rather the sum of all terms on the right hand sides of these two master equations 
vanish. This implies that population relaxation and thermal decoherence (in the case 
of finite temperature) are competing with the spin-spin interaction terms to create a 
steady state lcvcl of entanglement, as measured by the concurrence. 

We note, finally, that a somewhat different model system studied by S. 
Montangero, G. Benenti, and R. Fazio [3Ü| has found results for the pairwise 
concurrence that are somewhat similar to those we find for our system. Specifically, 
they have considercd the entanglement of a pair of spins within a qubit lattice in which 
there is disorder in the spin-spin couplings. They have identified a regime in which the 
pairwise concurrence is stable against such disorder in the couplings and has a value in 
the range of 0.2-0.3 We note that the numerical maximum for their "saturation 
value" of the concurrence is quite close to the analytical maximum we have derived 
in this paper (cf. Eq. I|28|l). It is interesting to observe that our analytic result for 
the maximum value of the steady state concurrence is constant for a range of system 
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parameters (cf. Eqs. (J32J-J2HJ)- Whether or not this analytical maximum holds also 
for other systems, such as the different one considered in |39|. is an open question. 

In summary, we have provided a detailed analytical and numerical analysis of 
decoherence for an interacting two-qubit model system having a Hamiltonian identical 
in form to that for the well-known two-qubit Heisenberg XY spin 1/2 system in 
the presence of an (effective) external uniform magnètic field. For T = 0, we 
have presented an analytic solution for the evolution of entanglement, measured by 
concurrence, for the case that both qubits are initially in their ground states; we 
have presented also numerical solutions for two other typical initial states. We find 
that our system is robust against decoherence: a steady state level of entanglement, 
controllable by the vàlues of the system parameters, is always reached for zero or finite, 
low temperatures. For the T = case, we have defined this steady state analytically 
and obtained the parameter vàlues that maximize its entanglement. For T > 0, the 
steady state level of entanglement is found to vanish at a finite temperature. Since 
our model interaction Hamiltonian describes also mesoscopic objects that interact via 
their spins (e.g., cf. Ref. HO]), it may be that a certain level of entanglement is robust 
against decohering interactions with an environment even for such objects. As noted 
by Ghosh et al. 01], even "the slightest degree of entanglement can have profound 
effects" on the properties of mesoscopic spin systems. 

We acknowledge stimulating discussions with Joseph H. Eberly, Hong Gao, Andrci 
Y. Istomin, Murray Holland, Ting Yu, and Peter Zoller. This work is supported in 
part by grants from the Nebraska Research Initiative and the W. M. Keck Foundation. 
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